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Shock  equation  of  state  of  multi-constituent  epoxy-metal  particulate 
composites 
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(Received  15  July  2010;  accepted  29  November  2010;  published  online  14  January  2011) 

The  shock  properties  of  epoxy-based  particulate  composites  have  been  extensively  studied  in  the 
literature.  Generally,  these  materials  only  have  a  single  particulate  phase;  typically  alumina.  This 
paper  presents  equation  of  state  experiments  conducted  on  five  epoxy-based  particulate  composites. 

The  shock  stress  and  shock  velocity  states  were  measured  for  five  different  composites:  two 
epoxy-aluminum  two-phase  composites,  with  various  amounts  of  aluminum,  and  three 
epoxy-aluminum- (metal)  composites,  where  the  metal  constituent  was  either  copper,  nickel,  or 
tungsten.  The  impact  velocities  ranged  from  300  to  960  m/s.  Numerical  simulations  of  the 
experiments  of  epoxy-Al  are  compared  with  mesoscale  simulations  of  epoxy-Al203  composites  to 
investigate  the  effect  of  the  soft  versus  hard  particulate;  additionally,  an  epoxy-Al-W  simulation 
was  conducted  to  investigate  the  material  properties  of  the  second  phase  on  shock  response  of  these 
materials.  In  these  epoxy-based  particulate  composites,  the  slope  of  the  shock  velocity-particle 
velocity  curve  appears  to  depend  on  the  epoxy  binder.  It  is  shown  that  the  addition  of  only  10  vol  % 
of  a  second,  denser  metallic  phase  significantly  affects  the  shock  response  in  these  composites. 

©  2011  American  Institute  of  Physics.  [doi:10.1063/l. 3531579] 


I.  INTRODUCTION 

Particulate  composites  comprised  of  a  polymer  binder 
with  one  or  multiple  particulate  phases  have  a  variety  of 
engineering  applications  with  their  high  strength  and  low 
density.  The  shock  properties  of  such  composites  are  of  in¬ 
terest  as  explosives  and  propellants  are  composed  of  particu¬ 
late  phases  in  polymer  binders.  There  are  several  studies  in 
the  literature  regarding  the  equation  of  state  of  single  particu¬ 
late  based  epoxy-Al203  (Refs.  1-5)  and  epoxy-WC  (Ref.  6) 
composites.  Although  these  single  component  systems  inter¬ 
act  in  a  complex  manner  with  shock  waves,  the  addition  of  a 
second  metal  or  ceramic  particulate,  such  as  in 
epoxy-Al -Fe203  (Refs.  7  and  8)  or  epoxy-Al-Mn02,^’*°  can 
result  in  even  more  complex  interactions.  The  propagated 
wave  observed  in  epoxy-Al203  broadens  at  low  input  stress 
due  to  the  increased  time  available  for  viscous  mechanisms.^ 
As  the  input  stress  increases,  the  epoxy-Al203  composite 
material  exhibits  viscoelastic  behavior.  Additionally,  the  re¬ 
lease  wave  velocity  is  a  strong  function  of  particle  velocity 
and  much  faster  than  the  initial  shock  wave.  ’  A  recent  in¬ 
vestigation  shows  that  epoxy-WC  composites  have  a  similar 
stress-strain  response  as  epoxy-Al203  composites,  as  well  as 
a  strong  dependence  of  release  wave  speed  on  particle 
velocity.® 

In  epoxy, **“*^  Carter  and  Marsh'*  observed  failure  of  the 
shock  velocity  to  extrapolate  to  the  ultrasonic  bulk  sound 
speed,  which  they  attributed  to  the  compressed  distance  be- 


“^Electronic  mail:  jennifer.jordan@eglin.af.mil. 


tween  the  polymer  chains  with  rigid  polymer  backbones.  Ad¬ 
ditionally,  they  observe  a  high  pressure  phase  transition,  at 
—23  GPa,  which  was  attributed  to  interchain  chemical 
reactions." 

This  paper  will  discuss  equation  of  state  experiments 
conducted  on  epoxy-based  particulate  composites.  These 
composites  contain  aluminum  powder,  for  two-phase  com¬ 
posites,  along  with  copper,  nickel,  or  tungsten  powder,  for 
three  phase  composites.  It  has  been  shown  that  the  matrix 
material  dominates  the  bulk  shock  response  of  epoxy-based 
particulate  composites  regardless  of  the  material  properties 
of  the  particulate  phase.”*  ®  The  third  phase  here  is  added  to 
investigate  the  effects  of  the  complementary  deformation  of 
the  two  powder  phases.  For  example,  it  was  shown  in  Ref. 
14  that  the  addition  of  a  hard  and/or  dense  third  phase 
changes  the  distribution  of  internal  energy  in  shocked  three 
phase  composites.  Stress  and  time  of  arrival  were  measured 
using  manganin  gauges.  The  effect  of  the  deformability  of 
the  particulate  is  investigated  through  simulations  of  the 
epoxy-Al  composite  studied  in  this  paper  with  the 
epoxy-Al203  composite  extensively  studied  in  the  literature. 
The  addition  of  a  second  particulate  phase  is  considered 
through  simulations  of  the  epoxy-Al-W  composite. 

II.  EXPERIMENTAL  PROCEDURE 

Five  materials  were  prepared  for  this  study — two  com¬ 
posites  containing  only  epoxy  and  aluminum,  at  two  volume 
fractions,  and  three  composites  containing  the  addition  of  a 
second  metallic  phase.  The  manufacturer  and  average  par- 
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TABLE  I.  Particulate  characteristics. 


Powder 

Supplier 

Average  particle  size 
(/rm) 

Aluminum  (X81) 

Toyal 

27 

Copper 

Atlantic  equipment  engineers 

37 

Nickel 

Atlantic  equipment  engineers 

44 

Tungsten^ 

H.  C.  Starck  (Kulite) 

37 

^Reference  15. 

tide  size  for  the  powders  are  given  in  Table  I.  The  appropri¬ 
ate  volume  fractions  of  powder  for  each  composite  were 
blended  into  Epon  826  and  cured  with  diethanolamine.  The 
composite  mixture  was  cast  into  blocks  and  nominally  34.9 
mm  diameter  by  3.2  mm  thick  samples  were  machined.  The 
density  of  each  composite  was  measured  using  pyncnometry. 
The  longitudinal  and  shear  ultrasonic  wave  speeds  in  the 
composites  were  measured  using  a  GE  Panametrics  25  HP 
Plus  ultrasonic  thickness  gauge  with  a  2.25  MHz  M106 
transducer  for  longitudinal  measurements  and  a  2.25  MHz 
VI 54  transducer  for  shear  measurements.  The  compositions 
of  all  the  materials,  along  with  their  measured  densities  and 
sound  speeds,  are  given  in  Table  II.  Micrographs  of  each  of 
the  composites  are  given  in  Pigs.  l(a)-l(e). 

Four  gas  gun-driven  equation  of  state  experiments,  with 
impact  velocities  ranging  from  300  to  960  m/s,  were  con¬ 
ducted  using  the  102  mm  diameter  single-stage  light-gas  gun 
at  NSWC-Indian  Head.  The  test  configuration^  used  in  this 
study  allows  for  the  simultaneous  measurement  of  three 
samples.  Each  sample  was  instrumented  with  two  50  fl 
manganin  gauge  (Dynasen,  Goleta,  CA),  which  have  an  ac¬ 
tive  gauge  area  of  40.3  mm^  and  a  nominal  thickness  of  0.05 
mm,  mounted  between  the  donor  plate  and  the  sample  and 
between  two  sample  disks.  An  aluminum  projectile  with  an 
aluminum  flyer  plate  (12.7  mm  nominal  thickness)  was  used 
to  impact  an  aluminum  driver  plate  (6.73  mm  nominal  thick¬ 
ness)  with  the  samples  mounted  on  the  rear  of  the  driver 
plate.  All  of  the  plates  were  machined  from  6061-T6  Al.  The 
gauge  mounted  between  the  driver  and  the  sample  provides 
and  “input”  stress  profile,  and  the  gauge  mounted  between 
the  two  sample  disks  provides  the  “transmitted”  stress  pro¬ 
file. 

III.  RESULTS  AND  DISCUSSION 

Three  impact  velocities  were  used  to  determine  the  re¬ 
sponse  of  these  epoxy-based  particulate  composites  to  shock 
loading.  epoxy-45Al  and  epoxy-Al-W  were  impacted  at  302 


FIG.  1.  Micrographs  of  (a)  epoxy-35  Al,  where  the  light  gray  particles  are 
aluminum;  (b)  epoxy-45Al,  where  the  light  gray  particles  are  aluminum;  (c) 
epoxy-Al-Cu,  where  the  light  gray  particles  are  copper  and  the  dark  gray 
particles  are  aluminum;  (d)  epoxy-Al-Ni,  where  the  light  gray  particles  are 
nickel  and  the  dark  gray  particles  are  nickel;  and  (e)  epoxy-Al-W,  where  the 
white  particles  ai'e  tungsten  and  the  dark  gray  particles  are  aluminum. 


m/s.  All  the  samples  were  impacted  at  the  identical  impact 
condition  of  620  m/s,  and  epoxy-45Al  was  also  impacted  at 
960  m/s. 

The  transmitted  pressure  pulses  for  each  sample  at  622- 
624  m/s  impact  velocity  are  shown  in  Fig.  2.  The  large  spikes 
in  the  epoxy-Al-Cu  and  epoxy-Al-Ni  traces  are  due  to  fail¬ 
ure  of  the  input  gauge.  The  regular  oscillations  on  several 
wave  profiles  appear  to  be  due  to  electronic  noise  rather  than 
oscillations  within  the  gauge  package.  There  is  little  observ¬ 
able  difference  in  the  wave  profiles  between  samples,  other 
than  the  magnitude,  which  may  be  explained  through  the 
measured  density  differences  of  the  materials.  The  wave  pro¬ 
files  show  rounding  at  the  peak  of  the  rise  indicating 


TABLE  II.  Material  compositions,  measured  densities,  and  measured  longitudinal  (Cl)  and  shear  (Cs)  sound 
speeds,  and  calculated  bulk  sound  speed  (Cj,). 


Material 

Density 

(g/cm^) 

Cl 

(km/s) 

Cs 

(km/s) 

Cb 

(km/s) 

Al 

(vol  %) 

Cu 

(vol  %) 

Ni 

(vol  %) 

W 

(vol  %) 

epoxy 
(vol  %) 

Epoxy-35Al 

1.725 

2.95 

1.49 

2.40 

35 

65 

Epoxy-45Al 

1.875 

3.13 

1.61 

2.52 

45 

55 

Epoxy-Al-Cu 

2.475 

2.77 

1.40 

2.25 

35 

10 

55 

Epoxy-Al-Ni 

2.513 

2.95 

1.49 

2.40 

35 

10 

55 

Epoxy-Al-W 

3.652 

2.39 

1.20 

1.95 

35 

10 

55 
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Time  (ns)  (km/s) 


FIG.  2.  (Color)  Transmitted  pressure  vs  time  traces  for  all  samples  at  622 
m/s  impact  velocity. 


FIG.  3.  (Color  online)  Shock  velocity  vs  particle  velocity  for  all 
experiments. 


dispersion.^  However,  the  wave  dispersion  does  not  appear  to 
depend  on  the  presence  or  type  of  the  second  phase  in  the 
material. 

In  each  experiment,  the  shock  velocity  (C/J,  shown  in 
Table  III,  was  determined  from  the  difference  between  the 
time  of  arrival  at  the  manganin  gauge  mounted  between  the 
donor  plate  and  the  sample  and  the  gauge  mounted  between 
the  two  samples  and  the  known  sample  thickness.  The  pres¬ 
sure  in  each  sample  was  determined  at  an  input  and  a  trans¬ 
mitted  location,  approximately  3.2  mm  apart,  by  averaging 
the  pressure  at  the  plateau  observed  in  the  manganin  gauge 
signals.  These  measured  pressures  are  also  presented  in  Table 
III  as  and  respectively.  In  several  experiments,  it 
appears  that  the  transmitted  pressure  is  higher  than  the  input 
pressure,  which  is  physically  impossible.  The  error  on  these 
average  measurements  is  approximately  10%,  which  equates 
to  the  input  and  propagated  pressures  being  equal.  The  par¬ 
ticle  velocity,  Up,  was  determined  using  the  measured  shock 
velocities  and  input  pressure  for  each  experiment  by  the  con¬ 
servation  equation 

P  =  PoU,Up,  (1) 

where  po  is  the  initial  density.  For  experiments  where  the 
input  pressure  was  not  measured,  the  transmitted  pressure 
was  used,  which  is  valid  based  on  the  experiments  where 
both  pressures  were  measured  showing  a  steady  wave  over 
the  3.2  mm  propagation  distance. 


The  shock  velocity  versus  particle  velocity  for  each  of 
the  samples  is  given  in  Fig.  3.  An  empirical  linear  equation 
of  state  can  be  fit  in  the  U^-Up  space 

JJ ^  —  Cq  -f  Slip,  (2) 

where  Cq  is  the  bulk  sound  speed  in  the  material  at  zero 
pressure  and  s  is  an  empirical  parameter.  For  the  epoxy-45  A1 
sample,  the  shock  velocity  was  measured  at  three  particle 
velocities,  in  addition  to  the  bulk  sound  speed  measured  us¬ 
ing  ultrasound,  and  a  least-squares  linear  fit  determined 

U,  =  2.52+  1.89up.  (3) 

Although  only  two  points  were  measured  for  epoxy- Al-W,  a 
similar  linear  fit,  with  the  measured  ultrasonic  bulk  sound 
speed,  was  constructed  resulting  in 

U,=  1.92  +  2.11up,  (4) 

which  may  be  compared  to  the  fit  of  the  epoxy-WC  compos¬ 
ite  with  46  vol  %  WC  from  Ref.  6  which  was  t/^=1.65 
+  2.\5Up.  It  is  interesting  that  epoxy-45 A1  and  epoxy- Al-W 
have  similar  behavior  to  composites  in  the  literature,"^^^  as 
seen  in  Fig.  4.  The  shock  velocity-particle  velocity  lines  are 
parallel  within  approximately  10%  of  each  other.  This  indi¬ 
cates  that  the  epoxy  binder  is  driving  the  shock  response  of 
these  materials,  with  small  variations  in  the  epoxy  resins  and 
curing  agents  possibly  accounting  for  the  10%  difference. 


TABLE  III.  Experimental  results  for  gas  gun  experiments  on  epoxy-based  particulate  composites,  where  is 
the  shock  velocity,  Pj,,  is  the  pressure  from  the  input  manganin  gauge,  is  the  pressure  from  the  transmitted 
manganin  gauge,  and  Up  is  the  calculated  paiticle  velocity. 


Sample 

Impact  velocity 
(m/s) 

u, 

(km/s) 

Pin 

(GPa) 

Pout 

(GPa) 

% 

(km/s) 

Epoxy-35Al 

622  ±8 

3.22 

3.3 

3.1 

0.593 

Epoxy-45Al 

302  ±2 

2.93 

1.05 

1.04 

0.190 

Epoxy-45Al 

622  ±8 

3.38 

3.4 

0.537 

Epoxy-45Al 

960  ±2 

3.93 

5.1 

0.693 

Epoxy-Al-Cu 

624  ±0.5 

3.03 

3.4 

3.6 

0.488 

Epoxy- Al-Ni 

624  ±0.5 

3.17 

3.8 

0.477 

Epoxy-Al-W 

302  ±2 

2.21 

1.2 

1.3 

0.153 

Epoxy-Al-W 

624  ±0.5 

2.82 

4.3 

4.4 

0.418 
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FIG.  4.  (Color  online)  Shock  velocity  vs  particle  velocity  for  all  experi¬ 
ments  compared  with  literature  data  for  epoxy-43  vol  %  AI2O3  (Refs.  2-5) 
and  epoxy-WC  (Ref.  6). 


For  the  materials  presented  here,  the  total  volume  frac¬ 
tion  of  particulate  is  the  same  in  epoxy-45Al  and  the  epoxy- 
Al-(metal)  composites,  with  epoxy-35Al  containing  the  same 
volume  fraction  of  aluminum  as  the  epoxy- Al- (metal) 
samples  and  10  vol  %  less  aluminum  than  epoxy-45  Al.  For 
epoxy-45Al  and  epoxy-Al-W,  the  ultrasonically  measured 
sound  speed  appears  to  he  linear  with  the  measured  Ug-Up 
data  points.  In  polymers,  the  extrapolated  Cq  value,  from  a 
linear  fit  to  the  experimental  shock-loading  points,  is  typi¬ 
cally  higher  than  the  ultrasonically  measured  bulk  sound 
speed. However,  Millett  et  al^  observed  that  this  behavior 
shifted  to  the  more  common  behavior  of  the  extrapolated  Cq 
correlating  with  the  ultrasonically  measured  bulk  sound 
speed  with  increasing  amounts  of  alumina,  which  is  observed 
in  these  materials. 

Millettt  et  al^  also  observed  that  the  slope,  s,  appears  to 
be  controlled  by  the  epoxy  binder,  since  there  was  little  dif¬ 
ference  in  the  slopes  of  two  alumina-epoxy  composites  with 
different  volume  fractions  of  alumina.  Despite  the  limited 
data  in  this  study,  this  also  appears  to  be  true  for  these 
epoxy-based  composites.  The  s  values  in  Eqs.  (1)  and  (2)  are 
higher  than  those  reported  by  Carter  and  Marsh  (1.51)," 
Munson  and  May  (1.66),"  and  Millett  et  al.  (1.47)  (Ref.  13) 
for  the  pure  epoxy  resin.  The  small  differences  between  the 


reported  literature  values  are  likely  due  to  differences  in 
hardeners  and  plasticizers  used  to  cure  the  epoxy  resin.  Ad- 
ditionally,  Setchell  and  Anderson  report  5=2.00  for 
43  vol  %  AI2O3  and  Millett  et  al.‘^  report  values  of  S 
=  1.63  for  their  “fully  loaded”  alumina-epoxy  composite  and 
s=  1.66  for  the  “half  loaded”  composite.  An  average  value  of 
s  fit  to  all  the  epoxy-43  vol  %  AI2O3  data^”^  is  1.84,  as 
shown  in  Fig.  3.  The  slopes  of  all  the  epoxy-based  compos¬ 
ites,  with  the  epoxy-43  vol  %  AI2O3  composite  showing  a 
shallower  slope  than  epoxy-45Al,  epoxy-Al-W,  and  epoxy- 
WC. 


IV.  NUMERICAL  MODELING 

Numerical  simulations  of  the  shock  response  of 
epoxy-Al  composites,  considered  in  this  study,  are  presented 
and  contrasted  with  simulations  using  epoxy-Al203,  which 
has  been  studied  extensively  in  the  literature,  to  investigate 
the  effects  of  material  properties  and  particle  morphology  of 
the  powder  within  the  matrix.  Results  of  three  component 
systems  using  epoxy-Al-W  for  comparison  with  Eq.  (4)  and 
to  investigate  the  effects  of  adding  only  10%  of  third,  dense 
material  in  numerical  simulations  are  also  presented. 

The  two-dimensional  microstructural  images  [similar  to 
Fig.  1(b)]  are  converted  into  vectorized  lists  of  particle 
boundaries  for  the  simulated  impact  of  the  domain.  Many 
investigations  use  randomly  distributed  circles,  squares  and 
cubes  with  appropriate  size  and  distribution  to  represent 
powder  particles  in  mesoscale  simulations®’"^"  with  excel¬ 
lent  results.  The  fine  details  of  the  particle  surface  are  usually 
not  as  important  as  nominal  aspect  ratio  (e.g.,  nominally 
round  or  flat")  and  the  size  distribution  due  to  the  fact  that 
small  asperities  on  the  particle  surfaces  may  be  diminished  in 
the  advection  step  in  hydrocodes.  The  use  of  digitized  micro¬ 
graphs  does  indicate  the  nominal  shape  and  size  distribution, 
but  the  intended  volume  (area)  fraction  of  constituents  must 
be  verified  to  account  for  the  fact  that  micrographs  are  a 
two-dimensional  representation  of  a  three-dimensional  mate¬ 
rial.  Thus,  the  resulting  digitization  is  highly  dependent  on 
the  precise  location  of  the  polished  reference  surface. 


s 


> 


FIG.  5.  (Color  online)  Initial  vectorized  microstructures  for  simulation  of  experiments  with  epoxy-45Al. 
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TABLE  IV.  Calculated  volume  fraction  and  density  of  digitized  microstractures  for  numerical  simulations  for 
epoxy-45  A1  and  epoxy- Al-W. 


Material 

Designation 

Volume 

fraction 

epoxy 

Volume 

fraction  A1 

Volume 

fraction  W 

Po 

(g/cc) 

Epoxy-45Al 

Microstructure  A 

52.8 

47.2 

0 

1.876 

Epoxy-45Al 

Microstructure  B 

51.8 

48.2 

0 

1.892 

Epoxy-45Al 

Microstructure  C 

53.1 

46.9 

0 

1.872 

Epoxy-Al-W 

Microstructure  D 

58.3 

29.1 

12.6 

3.877 

A  custom  MATLAB  program  was  used  to  threshold  the 
grayscale  values  of  the  particles,  separate  the  phases  within 
the  composite  and  locate  the  boundaries.  The  outline  of  the 
particles  was  then  converted  to  a  vector  list  representation 
and  scaled  from  a  pixel  to  length  scale  (similar  to  Ref.  18). 
The  program  created  an  input  file  for  each  specific  simula¬ 
tion  for  a  parallel  version  of  CTH,  which  is  a  finite  volume 
three-dimensional  shock  wave  physics  code.*® 

The  two-dimensional  microstructures  used  in  the  shock¬ 
loading  simulations  are  shown  in  Figs.  5(a)-5(c),  7,  and  8(a). 
The  microstractures  were  discretized  into  0.36  /xm  cells  in 
Figs.  5  and  8  and  0.06  /xm  cells  in  Fig.  7.  The  top  and 
bottom  boundaries  are  given  transmitting  boundary  condi¬ 
tions.  The  sides  of  the  computational  domain  are  prescribed 
as  symmetry  planes.  The  calculated  volume  fractions  and 
densities  are  given  in  Table  IV.  The  three  digitized  micro¬ 
structures  in  Fig.  5  were  each  impacted  at  302,  622,  and  960 
m/s  (as  shown  in  Table  III)  by  a  6061-T6  A1  half-space 
placed  above  the  sample  to  produce  an  infinitely  long  shock. 

The  sizes  of  the  micrographs  were  selected  based  upon 
simulations  (not  shown)  where  multiple  microstructures 
were  stitched  together  to  lengthen  the  sample  in  the  direction 


of  the  shock.  The  measured  difference  between  the  shock 
and  particle  velocity  was  within  the  difference  (<1%)  ob¬ 
served  by  changing  the  location  of  the  line  of  tracers  and 
using  different  digitized  microstractures. 

To  measure  the  shock  velocity  and  particle  velocity  for 
comparison  with  experiments,  a  horizontal  line  of  150  tracer 
particles  is  placed  within  the  composite  at  y= 0.002  cm  and 
y=0.022  cm  similar  to  the  methods  in  Refs.  6  and  20-22. 
The  shock  speed  was  calculated  by  determining  the  time 
when  the  pressure  of  each  tracer  particle  reached  half  of  its 
maximum  height.  The  shock  speed  is  then  calculated  using 
the  position  and  the  “half-maximum”  of  pressure  times  of  the 
tracers  that  were  initially  above  and  below  one  another.  All 
of  the  individual  speeds  are  then  averaged  together  resulting 
in  a  single  shock  speed  for  the  simulation.  The  particle  speed 
is  determined  by  taking  the  average  particle  speed  in  the 
direction  of  the  shock  of  the  tracer  particles  after  the  shock 
has  passed. 

The  results  from  these  simulations  are  shown  in  Fig.  6 
where  the  shock  velocity  is  shown  as  a  function  of  particle 
velocity  comparison  for  (a)  epoxy-45Al  and  (b)  epoxy-Al203 


[km/s]  [km/s] 


FIG.  6.  (Color)  Shock  velocity  vs  particle  velocity  comparison  for  (a)  epoxy-45Al  and  (b)  epoxy-Al203  composites,  (a)  The  black  squares  are  the  data  points 
from  Table  III  for  epoxy-45Al.  The  (2)  black  line  is  the  least-squares  fit  of  this  data:  {/j=2.45-t2.18Mp.  There  are  three  sets  of  blue  squares  that  represent  the 
results  from  taking  each  particle  morphology  (shown  in  Fig.  5)  and  subjecting  it  to  impacts  at  302,  622,  and  960  m/s.  The  (1)  blue  line  is  a  least-squares  fit: 
U,=2J6+2.03Up.  The  blue  “-t”  denotes  an  impact  of  the  microstructure  (a)  using  the  Hugoniot  data  from  Millett,  Bourne,  and  Deas  (Ref.  4).  The  blue  circle 
denotes  a  simulation  where  the  material  properties  are  the  same  as  for  the  blue  squares,  but  a  boundary  layer  interface  was  prescribed  between  the  powder 
particles  and  the  epoxy  matrix,  (b)  The  black  squares  are  taken  from  Ref.  4.  The  blue  squares  are  computed  similarly  to  those  in  (a)  but  the  results  are  taken 
from  a  digitized  micrograph  from  Fig.  1  in  Ref.  5  using  the  properties  of  epoxy  and  AI2O3.  The  blue  line  [denoted  by  (4)]  shows  the  least-squares  fit  to  this 
data:  U, =2.54+2. l2Up.  the  three  blue  “X”  points  show  the  results  using  the  parameters  from  Ref.  6.  The  three  points  are  taken  from  a  “stitched”  micrograph 
of  epoxy-45 A1  where  the  A1  particles  are  given  the  material  properties  of  AI2O3. 
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TABLE  V.  Parameters  for  the  material  equation  of  state. 


Material 

EOS 

Po 

(g/cc) 

<^0 

(km/s) 

■«i 

V 

Epon-828 

Mie-Griineisen 

1.14 

2.69 

1.51 

0.35 

Al 

Tabular 

2.71 

5.25 

1.37 

0.28 

AI2O3 

Mie-Griineisen 

3.98 

7.93 

1.20 

0.28 

Al  6061-T6 

Mie-Griineisen 

2.71 

5.22 

1.37 

0.28 

W 

Mie-Giiineisen 

19.24 

3.98 

1.24 

0.28 

composites.  The  parameters  used  for  the  equation  of  state 
and  strength  models  are  shown  in  Tables  V  and  VI. 

In  Fig.  6(a),  the  three  black  squares  are  the  experimental 
data  points  from  Table  III  for  epoxy-45Al.  The  black  line 
labeled  (2)  is  the  least-squares  ht  of  this  data:  f/j=2.45 
-i-2.18Mp.  There  are  three  sets  of  three  blue  squares  above  the 
experimental  line  that  represent  the  results  from  taking  each 
particular  particle  morphology  (shown  in  Fig.  4)  and  subject¬ 
ing  it  to  impacts  at  302,  622,  and  960  m/s.  The  simulations 
are  repeated  for  each  impact  condition  with  different  micro¬ 
structures.  The  blue  line  labeled  (1)  is  a  least-squares  fit  of 
the  numerical  data  points:  Us=2.16  +  2.03up.  The  blue  “-I-” 
denotes  an  impact  of  the  microstructure  “A”  using  the  epoxy 
Hugoniot  data  from.^  Note  that  changing  the  Cq  and  values 
for  the  epoxy  only  result  in  an  averaged  decrease  in  the 
shock  speed  (and  slope  which  is  not  shown)  for  a  particle 
velocity  of  302  m/s.  The  two  red  diamonds  are  taken  from 
the  epoxy-Al-W  data  in  Table  III.  The  red  line  denoted  by 
curve  (3)  is  Eq.  (4).  The  red  square  above  this  line  is  the 
simulated  result  using  the  microstructure  from  Fig.  8. 

Figure  6(b)  presents  the  shock  speed  (11^)  and  particle 
speed  (up)  data  from  the  fully  loaded  AI2O3  samples  in  Ref. 
4  (black  squares)  and  Ref.  3  (black  “X”  data).  The  blue 
squares  are  computed  similarly  to  those  in  Fig.  5(a),  but  the 
results  are  taken  from  simulations  using  the  particle  mor¬ 
phology  shown  in  Fig.  1  in  Ref.  5,  shown  here  in  Fig.  7.  The 
blue  line  [denoted  by  (4)]  shows  the  least-squares  fit  to  this 
data:  Us  =  2.54+2.l2up.  The  three  blue  X  points  show  the 
results  using  the  epoxy  parameters  from  Ref.  6.  The  three 
points  are  taken  from  a  “stitched”  micrograph  of  epoxy-45Al 
where  the  A1  particles  are  given  the  material  properties  of 
AI2O3.  The  “stitched”  micrograph  was  twice  as  long  as  those 
shown  in  Fig.  4  and  the  three  points  are  U^,  Up  pairs  mea¬ 
sured  at  three  different  locations  for  lines  of  tracer  particles 
placed  at  four  different  locations.  The  U^,  Up  pairs  differ  by 
about  1%  in  shock  speed  and  in  particle  speed. 

The  results  shown  in  Fig.  6(a)  indicate  that  changing  the 
particle  morphology  (e.g.,  see  Fig.  5)  does  not  result  in  an 
appreciable  change  in  the  U^,  Up  values.  This  is  also  true  for 


simulations  of  epoxy-Al203  samples,  as  seen  in  Fig.  7,  where 
the  morphologies  are  drastically  different  not  only  in  particle 
distribution,  but  size  and  shape.  It  is  interesting  that  the 
slopes  of  the  epoxy-45AI  simulations  match  very  well  be¬ 
tween  simulations  and  experiments  and  that  the  experimental 
data  for  the  epoxy-Al203  is  between  literature  data  depend¬ 
ing  on  the  material  parameters  used.  It  is  most  probable  that 
uncertainties  in  the  density  and  amount  of  fill  of  the 
epoxy-Al203  sample  were  the  main  cause  of  the  discrepancy 
between  the  experiments  and  simulations.  This  is  supported 
by  the  fact  that,  regardless  of  the  type  of  material  that  is  put 
into  the  epoxy  (Al,  AI2O3,  etc.),  it  has  been  suggested  that 
the  epoxy  matrix  drives  the  response  of  the  materials  when 
the  volume  fraction  of  the  hll  material  is  low  enough  that 
particle-particle  interactions  do  not  have  a  strong  effect.^’® 
The  two  red  diamonds  in  Fig.  6(a)  are  taken  from  the  epoxy- 
Al-W  data  in  Table  III.  The  red  line  denoted  curve  (3)  is  Eq. 
(4).  The  red  square  above  this  line  is  a  simulated  result  using 
the  microstructure  from  Eig.  8  impacted  at  450  m/s  resulting 
in  values  of  Up  =  332  m/s,  I7j=2.921  km/s.  The  image  from 
Fig.  1(e)  was  rotated  90°  for  the  simulation  in  Fig.  8  to 
elongate  the  material  in  the  direction  of  shock.  The  pressure 
distribution  is  shown  in  Fig.  8(b)  at  the  stage  where  the 
shock  has  almost  reached  the  bottom  of  the  material  and  it  is 
clear  that  the  shock-front  is  nonuniform  (the  colors  corre¬ 
spond  to  blue=0  GPa,  and  red =2  GPa).  The  W  particles 
clearly  disrupt  the  planarity  of  the  shock  front,  but  do  not 
affect  the  slope  of  the  U^,  Up  curve  as  indicated  in  Fig.  6(a) 
and  from  Eq.  (4).  The  W  particles  are  much  less  deformable 
than  epoxy  or  the  surrounding  Al  particles,  which  undergo 
much  more  bulk  deformation  along  with  the  matrix. 

V.  CONCLUSIONS 

This  paper  discussed  equation  of  state  experiments  and 
complementary  mesoscale  simulations  conducted  on  five 
epoxy-based  particulate  composites.  Two  of  the  composites 
contained  only  aluminum  (35  and  45  vol  %)  and  epoxy.  The 
remaining  three  composites  contained  35  vol  %  aluminum 


TABLE  VI.  Parameters  for  material  strength  models.  Al  and  AI2O3  used  Johnson-Cook  models  built-in  to  CTH 
material  libraries. 


Material 

Strength  model 

A 

(MPa) 

B 

(MPa) 

C 

m 

n 

T,„ 

(K) 

Epon-828 

Johnson-Cook 

84 

15 

0.01 

0.53 

0.20 

350 

AI2O3 

Johnson-Cook 

3.98 

7.93 

1.20 

0.28 

0.60 

2073 

w 

Johnson-Cook 

1510 

177 

0.016 

1 

0.12 

1748 
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Materials  at  O.OOe+00  seconds 


X  (urn) 

FIG.  7.  (Color  online)  Digitized  migrograph  taken  from  Ref.  5.  The  com¬ 
puted  density  is  2.496  g/cm^  which  corresponds  to  48  vol  %  of  epoxy  and 
52%  AI2O3. 


with  10  vol  %  of  an  additional  metal  (copper,  nickel,  or 
tungsten)  in  an  epoxy  binder.  The  pressure  and  shock  veloc¬ 
ity  of  these  epoxy-based  composites  were  measured  at  im¬ 
pact  velocities  ranging  from  300  to  960  m/s,  using  manganin 
gauges  for  stress  and  time  of  arrival  measurements.  Numeri¬ 
cal  simulations  have  been  conducted  comparing  the 
epoxy-Al  from  this  study  with  epoxy-Al203  previously  pub¬ 
lished  in  the  literature  to  investigate  the  effect  of  the  material 
properties  of  the  particulate.  Additionally,  a  mesoscale  simu¬ 
lation  of  epoxy-Al-W  has  been  conducted  to  determine  the 
effect  of  a  second  particulate  phase. 

The  epoxy  binder  was  found  to  drive  the  behavior  of 
these  materials.  The  experimental  pressure  versus  time  traces 
as  measured  with  manganin  pressure  gauges  showed  little 
difference  between  the  hve  materials  studied.  Comparing  the 
shock  velocity  versus  particle  velocity  curves  for  the  five 
materials  in  this  study  as  well  as  two  from  the  literature 
revealed  nominally  linear  behavior  over  the  range  studied 
and  parallel  linear  hts  varying  only  with  sample  density.  Me¬ 
soscale  numerical  simulations  of  comparable  formulations 
containing  aluminum  and  alumina  particulates  demonstrated 
that  the  materials  properties  of  the  particulate  did  not  greatly 
affect  the  shock  behavior,  excluding  density  effects.  Simula¬ 
tions  of  epoxy-Al-W  showed  that  the  tungsten  greatly  af¬ 
fected  the  shock  fron  at  the  mesoscale,  but  this  behavior  is 
not  manifested  in  the  bulk  shock  behavior. 
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FIG.  8.  (Color)  Digitized  micrograph  of  Fig.  1(e).  The  image  was  rotated 
90°  to  elongate  the  material  in  the  direction  of  shock.  [t/iinpact=450  m/s, 
w^=332  m/s,  Us=2.921  km/s,  which  is  shown  in  Fig.  6(a)].  The  pressure 
distribution  is  shown  at  the  stage  where  the  shock  has  almost  reached  the 
bottom  of  the  material.  The  shock-front  is  nonuniform.  The  colors  corre¬ 
spond  to  blue=0  GPa,  and  red=2  GPa.  The  temperature  distribution  is 
shown  in  the  bottom  part  of  the  hgure.  The  W  particles  are  clearly  distin¬ 
guished  from  the  soft  A1  particles,  which  undergo  much  more  bulk  defor¬ 
mation  along  with  the  matrix. 
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